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ABSTRACT 

In the upcoming synoptic all-sky survey era of astronomy, thousands of new multiply imaged quasars 
are expected to be discovered and monitored regularly. Light curves from the images of gravitationally 
lensed quasars are further affected by superimposed variability due to microlensing. In order to 
disentangle the microlensing from the intrinsic variability of the light curves, the time delays between 
the multiple images have to be accurately measured. The resulting microlensing light curves can then 
be analyzed to reveal information about the background source, such as the size of the quasar accretion 
disc. In this paper we present the most extensive and coherent collection of simulated microlensing 
light curves; we have generated > 2.5 billion light curves using the GERLUMPH high resolution 
microlensing magnification maps. Our simulations can be used to: train algorithms to measure lensed 
quasar time delays, plan future monitoring campaigns, and study light curve properties throughout 
parameter space. Our data are openly available to the community and are complemented by online 
eResearch tools, located at http://gerlumph.swin.edu.au 

Subject headings: gravitational lensing: micro - accretion, accretion discs - quasars: general - astro¬ 
nomical databases: miscellaneous - virtual observatory tools 


1. INTRODUCTION 

Gravitationally lensed quasars are unique natural lab¬ 
oratories for exploring the physics of extreme environ¬ 
ments in the Universe. Microlensed quasars in particular 
represent the only systems known so far where properties 
such as the temperature profile of the quasar accretion 
disc can be probed. 

Quasar microlensing was first observed as uncorrelated 
light curve variability attributed to compact stellar mass 
objects within the galaxy-lens that lie close to the line 
of sight (llrwin et al. 119891: iGhang fc Refsdallll979l. 11984 


Paczvnskil 119861 : iKavser et al .1 1 198 (it ISchneider fc Weiss! 
1987fl . From this it has been possible to derive con¬ 


straints on the accretion disc size and temperature profile 
in a number of systems. 

Accretion disc size estimates for wavelengths be¬ 
tween 0.07 and 2 microns vary in the range 14.5 < 
log(ri/ 2 ) < 17.3, where ri/ 2 is the half light ra- 
dius measured in cm (e.g. see fPoolev et, all 120071: 


Morgan et al Il20l(l iBlackburne et al.ll2011l : lMunoz et al.l 


20111 : iJimenez-Viccnte et al.ll2014L and references therein 

for results on 25 systems). Temperature profiles are also 
found to be rather flat, i.e . independent of wavelength 
(e.g. see iBlackburne et al.l 120111 : 1.Iimenez-Vicente et all 
1201411 . Thes e results are in disagreeme nt with the thin 
disc theory (jShakura fc Sunvaevl ll973 ). which predicts 
sizes down to an order of magnitude smaller and tem¬ 
perature profiles with a 4/3 power law dependence on 
wavelength. 

Another experiment that can be performed using 
gravitationally lensed quas ars is measuri ng the value 
of Hubble’s constant, Hq f R.cfsdall 1964j). This has 
been perf ormed in at least 2 4 multiply imaged sys¬ 
tems (see iTortora et all 12004 iParaficz fe Hiorthl 120101 
lEulaers fe Magainll2011l : iRathna Kumar et alJl2014L and 

references therein) to yield values in the range 50 — 100 


km s _1 Mpc -1 . The wide range of values, which is less 
accurate than results obtained from ot h er methods (e.g. 
Ceph eid variability; iRiess et ahl 120111 : iFreedman et all 
120131 . is d ue to uncertainties and degeneracies of the lens 
potential (lKocha nekll2003 iKochanek fc Schechte rl 12004 

IQguril[2007t ISchneider fc Slusell2013HSuvu et al.ll2014fl . 

In order to obtain light curves suitable for Hq and 
accretion disc studies, lens systems have to be moni¬ 
tored for periods of months, or even years (e.g . see 
Gotti 1 198 it iKundic et al.l 119971 : lEigenbrod et al.l 120051 : 
Fohlmeister et al.112008H . Moreover, accurate observa¬ 
tional constraints, such as the positions and photom¬ 
etry of the multiple images, are required to produce 
accurate models of the galaxy-lens. Extracting the 
correct time delays can be further complicated by the 
onset of microlensing, which introduces additional un- 
cor related variability between the multip le images (e.g. 
see iHoiiati et al.l 120131 : iTewes et al1l2013lh In fact, mi¬ 
crolensing is expected to be taking place in all gravita¬ 
tionally lensed quasars, although with varying strength. 

The effect of microlensing is usually modeled using a 
magnification map: a pixellated map of the magnifica¬ 
tion in the source plane induced by the foreground mi¬ 
crolenses (see Figure [lj. Such maps are produced us¬ 


ing different im plementations of the i nverse ray-shooting 
technique (e.g. IKavser et all 11983 IWambsganssI H9 99: 


Kochanekll2064l: Thompson et alJl201Ct iMediavilla et al l 

20111 : 1 Vernardos fc Flukell2014alh 


Due to the relative motion of observer, lens, and 
source, the source appears to move across the map 
crossing regions of high and low magnification and pro¬ 
ducing a variable light curve. Simulated light curves 
have been used to study high magnification e vents 
(1 Wittl Il990t lYoneharall20dH iShalvapin et al.l 120021). au¬ 
tocorrelation functions ([Seitz et al.111994 iLewis fc Irwinl 
ll996HWvithe fc Loebll2002lh and derivative distributions 
( Wvithe fc Websteilll999H . iKochanekl ()2004f l developed a 
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quantitative Bayesian analysis technique which h as been 
used t o study the back ground quasar (e.g. iMorgan et al.1 

a I Da i et al.1 1 2 0 1 til and the foreground galaxy-lens 
Mosaucra ct al.ll2013D . 

From the ~ 90 known multiply imaged quasars 
(iMosaucra fc Kochanekl 12011 ). ~ 25 have been stud¬ 
ied using microlensing techniques, either as single ob¬ 
jects or in groups of a few. However, this is about to 
change due to the upcoming synoptic all-sky survey fa¬ 
cilities, like the L arge Synoptic Survey Telescope (LSST; 
ILSST Science Collaboration et all 120091) , which are ex¬ 
pected to discover thousands o f new lensed quasars suit - 
ablc for microlensing studies (iQguri fe Marshall!12010ft . 
Moreover, the high cadence of observations from these 
instruments will provide nearly effortless monitoring. 
These data are exp ected to increase the accuracy of Hn 
meas urements Ce.g. lCoc fc Moustakasll2009t iDobkc ct all 
120091) , and improve our c urrent techniques for constrain¬ 
ing quasar structure (e.g. lSluse fe Tew es 2014). It is cru¬ 
cial for microlensing to move from single object to pa¬ 
rameter space studies. 

As a theoretical counterpart of future quasar 
microlensing observational campaigns, the Graphics 
Processing Unit-Enabled High Resolution cosmolog¬ 
ical MicroLensing param eter survejQ (GERLUMPH; 
IVernardos fe Fluka l2014al) has already generated > 
70,000 microlensing magnification maps, the largest and 
most complete collection yet produced. The parame¬ 
ter space of convergence k, shear 7 , and smooth mat¬ 
ter fraction s, (see next section) is covered in unprece¬ 
dented detail, allowing for co mprehensive explorations o f 
microlensing properties (e.g. IVernardos fe Flukdl2014bl) . 

In this paper we present how we used the GERLUMPH 
maps to generate > 2.5 billion simulated microlensing 
light curves, the largest and most extensive set of light 
curves available to date. Our approach spreads the mod¬ 
elling effort all over the parameter space rath er than fo¬ 
cusing it on a single object, e.g. the method of lKochanekl 
(12001 may extract up to 10 6 light curves from a single 
source profile and magnification map, while we are ex¬ 
tracting 2000 light curves per source profile from 51,127 
maps covering the parameter space (see next section). 
Hence, our data is not designed for detailed modelling 
of individual objects but it is designed for studying the 
robustness of, and degeneracies in, models for many in¬ 
dividual objects across the parameter space. Using these 
simulations to unveil systematic errors introduced in the 
modelling process will hopefully lead to better measure¬ 
ments of H 0 and accretion disc constraints. 

We present our approach to extracting simulated mi¬ 
crolensing light curves in Section [2] Our data are de¬ 
scribed in Section [3] and are openly accessible for down¬ 
load, complemented by online analysis tools that we in¬ 
troduce in Section |4] Finally, we conclude our paper and 
present future prospects in Section [5] 

2. APPROACH 

We have used 51,227 microlensing magnification maps 
from the G ERLUMPH online resource , generated using 
the GPU-D (|Vernardos fe Flukel 120144 iThompson et al.l 
l 20 ll direct Graphics Processing Unit (GPU) im¬ 
plementation of the inverse ray-shooting technique 

1 http://gerlumph.swin.edu.au 


(|Kavser et al.lll986h . The map resolution is 10,000 pixels 
on a side, the map width is set to 25 Einstein radii (i?Einj 
see below), and the mass of the microlenses are fixed at 
1 Mg. The maps are extracted from the region of pa¬ 
rameter space with 0 < k < 1 and 0 < 7 < 1.3, which 
contains most of the models of the ga laxy-lens of the 
currently known systems (see figure 2 of lVernardos et all 
12014 or see the relevant GERLUMPH online too0). 

The smooth matter fraction is defined as s = k s /k, 
where k s is the contribution to the total convergence 
by the smooth matter component. For each k, 7 com¬ 
bination we use 10 maps with different smooth matter 
content: 0 < s < 0.9 in steps of 0.1. The collection of 
maps used in this study is contained within the slightly 
l arger s et of maps used previously in IVernardos fc Fluka 
(l2014bft . Properties of the GERLUMP H hig h resolution 
maps are examined in lVernardos et ahl (12014ft , while a de¬ 
scription of the GER LUMPH data, infra s tructur e, and 
tools, can be found in lVernardos fc Flukd (:2014a). 

The characteristic microlensing scale length in the 
source plane is the Einstein radius: 


Apiii 


D os D ls 4 G(M) 
D 0 \ c 2 


( 1 ) 


where D 0 \, D os , and D\ s , are the angular diameter 
distances from observer to lens, observer to source, 
and lens to source respectively, (M) is the mean mass 
of point-mass microlenses, G is the gravitational con¬ 
stant, and c is the speed of light. A typical range of 
values for i?Ein can be obtain ed from the sample of 
87 lensed quasars compiled by iMosauera fc Kochanekl 
<2QHD: 5.11 ± 1.88 x 10 16 cm. Although these authors 
have used estimates of the lens redshift in a number of 
cases, their result is consistent with the CASTLE80 sam¬ 
ple of 59 systems with both lens and source re dshifts mea¬ 
sured ('5.35±2.20xl0 16 cm. iFalco et al.ll2001ft . I n the fo l¬ 
lowing , we use the mean from the IMosauera fc Kochanekl 
( 2 OII I sample as a typical value for i?Ein- This leads to 
a pixel size of the high-resolution GERLUMPH maps of 
~ 1.28 x 10 14 cm. 

Because the microlensing effect depends weakly on the 
shape of the underlying accretion disc brightness profile 
dMortonson et al.l 120051) . the half light radius, 77 / 2 , is a 
convenient way to parametrize a range of disc models 
consistently. We model the quasar source profile as a 
face-orQ Gaussian disc, i.e. / (r) = exp (— r 2 /2a 2 ) , of 
varying size a that is related to 77/ 2 through 

r \/ 2 = 1 . 18 < t . ( 2 ) 

We truncate the profile at r = 3(7, which is the radius 
containing 99.7% of the total brightness. Therefore, the 
total size, or diameter, of the profile, d, which is used 


2 http://gerlumph.swin.edu.au/macromodels/ 

3 http://www.cfa.harvard.edu/castles/ 

4 Quasars are thought to be viewed mostly face—on within the 
curr ent framework o f the unified model for Active Galactic Nu¬ 
clei (Antonucci 1993; Urry & Padovani 1995). This is supported 
by microlensing results for a t least one system (Q2237+0305; 
IPoindexter fc Kochanekl 120 1 Oi l- However, depending on the (un¬ 
known) width of the dusty torus surrounding the central regions 
of the quasar, the viewing angle could be very different to face on. 
Since microlensing behaviour depends only weakly on shape, a disc 
that is not viewed face-on will just appear smaller. 
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Figure 1 . Central 2000-pixel (5 i?Ein) wide region of a magni¬ 
fication map with (^, 7 , s) = (0.52,0.36,0.5), produced by 5797 
microlenses. The locations of the sampled light curves and their 
intersection points are shown in grey. The magnification along the 
thick (green online) curve is shown in Figure [2] for 25 different 
source profiles. 

to determine the size of the effective map (see below), is 
equal to 6cr. We consider 25 profiles of varying size d: 
from 2 x 10 15 cm to 2 x 10 16 cm (~ 0.8 to ~ 8 light days) 
in steps of 2 x 10 15 cm, and from 2 x 10 16 cm to 1.7 x 
10 17 cm (~ 8 to ~ 65 light days) in steps of 1 x 10 16 cm. 
Consequently, we cover the range 14.6 < log(r 1 / 2 ) < 
16.5 which contains most of the current microlensing size 
estimates of the accretion disc (see Section [TJ. 

As it will be shown in the next section, includ¬ 
ing wider profiles produces large amounts of additional 
data without significant change in the light curves. 
Such a behaviour is expected because the microlens- 
ing effect is more prominent for smaller sources (e.g. 
IWambsganss fe Paczvnskilll991h . while the magnification 
of large sources tends to the macro-magnification (see 
Equation [4]). 

Light curves extracted from the original GERLUMPH 
maps hold magnification information for “pixel-size” 
sources, i.e. sources that are smaller than the pixel size 
of the maps. To get information for a source profile with 
finite size one has to convolve it with the maps. This 
can be achieved using the convolution theorem and a 
GPU to accelerate compu ting Fourier transforms (see 
IVernardos fc Flukcl l2014ai for a detailed description of 
this technique). Assuming the map is periodic, which 
is required by the Fourier transform, leads to spurious 
magnification values around the edges of the convolved 
map. The size of these regions is equal to half the source 
profile size, i.e. d/2. The largest profile we use is 1332 
pixels wide, therefore, by disregarding areas of 700 pixels 
around the edges of the map we are not affected by spu¬ 
rious magnification values. The resulting effective con¬ 
volved maps, from which the light curves are extracted, 
have a width of 8600 pixels, or 21.5i?Ein- 

We extract 2,000 light curves from each effective map. 


Ideally, this number would be as high as possible, how¬ 
ever, increasing it further will produce much more data 
than what we can currently store. It turns out that this 
number is sufficient to adequately sample the underlying 
magnification probability distribution of the magnifica¬ 
tion maps (see next Section). 

The light curve locations are randomly selected once, 
and then used for all the maps. The microlens positions 
for generating the GERLUMPH maps are set randomly 
using different random seeds for each map. Therefore, 
the locations of caustic networks between maps are not 
correlated and the fixed light curve locations do not bias 
the actual light curves extracted. This approach reduces 
the number of parameters that we have to keep track of, 
and avoids issues with generating random numbers, e.g. 
in case of reproducing our results using different compiler 
versions. In Figure [U we show a 2000-pixel-wide region 
from the center of a fiducial magnification map with the 
locations of the extracted light curves. 

The light curve length is set to 1.5i?Ein and sampling 
length to 0.0025i?Ein (this is set by the resolution of 
our maps); using th e values fo r Ryjn and the effective 
source velocity from iKochancO (j2004f ) for Q2237+0305 
the light curve would be ~ 22 years long, sampled every 
~ 13 days. Ideally, longer light curves would be desired, 
however, further increasing the light curve length will 
produce much more data than what we can currently 
store. A length of 1.5i?Ein is adequate to capture the 
microlensing variations (i?Ein is the typical scale length 
for the onset of microlensing) and corresponds to ~ 3er of 
the Gaussian distribution we used to model our largest 
source profile. 

The actual number of pixels crossed varies according to 
the orientation of the light curve: 600 pixels are crossed 
in the horizontal or vertical direction, while only 425 pix¬ 
els are crossed at an angle of 45°. Obviously, in the latter 
case some pixels are sampled more than once leading to 
the appearance of short flat parts, or steps, in the light 
curve. This can also arise due to the orientation of the 
light curve with respect to a caustic and is a pathologi¬ 
cal behaviour due to the finite resolution of a pixellated 
magnification map. Some interpolation prescription be¬ 
tween neighbouring pixels would resolve this and make 
the light curves smoother. However, we currently avoid 
making any assumptions on such smoothing procedures 
and provide the raw sampled data. 

Our approach to pixel sampling lies between Bresen- 
hams’ algorit hm (lBrese nhamlll965l) and a supercover al¬ 
gorithm (e.g. I Andres et al.l 1 199711 . which both belong to 
the digital differential analyzer (DDA) class of algorithms 
used in computer graphics for rasterization of geometri¬ 
cal shapes. In fact, infinitely reducing the scaling length 
would lead exactly to a supercover pixel sampling, i.e. 
sampling of all the pixels a line crosses. The maxi¬ 
mum amount of light curve information extracted from a 
magnification map would be possible by calculating each 
portion of a curve that lies in a given pixel. However, 
our choice of sampling avoids the need for complex data 
structures, increased data size (two values required for 
each pixel), and resembles more the observational data 
which inevitably produce discrete samples. 

2.1. The data 
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We have carried out >1.28 million convolutions across 
5 days on a GPU supercomputer. Each combination of 
map, profile, and l?Ein is stored in an indexed database; 
the index of each entry points to a directory in a flat file 
system that holds all output (more details on the choices 
we made fo r the backen d database an d file storage can 
be found in IVernardos fc Flukdl2014ah . Each light curve 
consists of 600 sampled magnification values which are 
represented as 32-bit floats. The 2000 light curves per 
convolution are stored in a binary file of ~ 4.6 MB in 
size. We have generated a total of ~ 2.5 x 10 9 light 
curves, corresponding to 5.6 TB of data. All our data 
are freely accessible to analyze and download using the 
eTools described in Section 0J 

3. RESULTS 

As a first test, we consider maps for pixel-size sources 
and extract the magnification probability distribution 
(MPD) based purely on the light curve magnification val¬ 
ues from 2000 light curves. Then, we compare this light 
curve MPD to the full MPD extracted from the same 
map region using the Kolmogorov-Smirnov test with the 
null hypothesis that the two distributions are the same. 
We find that ~ 8% of the 51,227 tests fail with a p- 
value < 0.05. We note that the light curves intersect 
at 5,995 pixel locations (same for every convolved map), 
which are counted more than once in the calculation of 
the MPD, but given the low number of failed tests we 
choose to ignore this. 

More than 95% of the failed tests occur for maps pro¬ 
duced by < 10 4 microlenses. This means that the caus¬ 
tic networks are less dense and a higher number of light 
curves would be needed to better probe the underlying 
map MPD. Indeed, increasing the number of light curves 
to 4000 and 8000 leads to ~ 6.5% and ~ 4% of failed tests 
respectively (but increases the size of the data products 
by factors of 2 and 4). Finally, it is expected that when a 
pixel-size source map is convolved with a source profile, 
the caustic networks, and the corresponding MPD, are 
smoothed out. We perform the same MPD comparisons 
using 2000 light curves from convolved maps, and find 
that even the smallest source profile, i.e. d = 2 x 10 15 cm 
or ~ 0.04i?Ein, leads to < 1% of failed tests. 

The Kolmogorov-Smirnov test is only one method of 
assessing the usefulness of any sub-sample from a pop¬ 
ulation. The test is more sensitive to differences near 
the peak of the distribution, rather than at the extremes 
where the i ndividu al probabiliti es are low (see section 
4.1 of IVernardos fe Flukd I2014bi ). We reiterate that it 
is straightforward to generate additional lightcurves per 
map, and indeed more independent magnification maps, 
but it is non-trivial to store them. 

In Figure [2] we plot the magnification values along 
the trajectory shown in Figure [l] for 25 different pro¬ 
file sizes. We convert the magnification of the unlensed 
source luminosity into magnitude change with respect 
to the macro-magnification produced by the galaxy-lens 
without the additional effect of microlensing: 

( T total \ / t macro 

—— ) — 2.5 log (—-— 

= 2.5 log (n) - 2.5 log (/x t h) 

= 2.5 log \ , (3) 

\ Mth / 



Figure 2. Magnification values, expressed in units of Amag 
(Equation [3}, along the trajectory shown in Figure^] The thick 
black line (magenta online) corresponds to the smallest source 
size for which the microlensing induced fluctuations are the most 
prominent. The double peaks disappear for a profile with a size 
of 0.4i?Ei n , roughly equal to the double peak separation. The 
pixel-size source light curve is shown in grey. 


where L is the unlensed source luminosity, L total is 
the total magnified flux, L macro is the magnified flux 
due to the galaxy-lens only, and the theoretical macro¬ 
magnification is: 


1 

Mth ~7Z To o ■ 

(1 — kY — 


(4) 


For small source sizes, the light curve has a characteris¬ 
tic double-peaked shape, with the peaks corresponding 
to the caustic crossing events shown in Figure [Q As 
the profile size gets larger the two peaks are smoothed 
out and disappear for the profile with size 2 x 10 16 cm, 
or ~ 0.4i?Ein, which corresponds roughly to the length 
of the light curve that lies within the caustic. For even 
larger profiles the light curve becomes almost flat and 
microlensing fluctuations are much less prominent as ex¬ 
pected. 

In the original implementation of GPU-D 
dThompson et al.l l2010f) . we used a random distri¬ 
bution of light rays. Compared to the alternative, i.e. a 
regular grid, this was a compromise that helped amortize 
the computational cost of additional GPU-computation 
while allowing us to explore billion lens configurations 
with abitrary numbers of source-plane pixels. To this 
end, light ray positions were generated on the CPU 
(host memory) while the light ray deflections were 
being calculated on the GPU. As a result, there is an 
uncertaininty in the magnification values from each map 
pixel: counting IVjj rays in a given pixel is accompanied 
by an approximately yJN^j error, i.e. Poisson-like. The 
code comparisons made by iBate et abl (201 (1) showed 
that the these per-pixel magnification errors were small, 
and become less significant as the maps are convolved 
with realistic source profiles. 

Low magnifications correspond to low ray counts and 
consequently larger errors with respect to the value itself 
(i.e. V -1 / 2 ). Therefore, we expect larger fluctuations in 
the low magnification parts of a simulated light curve for 
sources smaller than the size of the map pixels. Indeed, 
such a behaviour is observed for the pixel-size source 
light curve shown in Figure [2] As soon as the pixel 
size source map is convolved with a profile - even the 




















2.5 billion microlensing light curves 


5 



Amag 

Figure 3. Magnification errors in units of Amag (same as y-axis 
of Figure |2j. The exact error for convolved magnification values 
(Equation [ 6 j is computed in 100 locations within a magnification 
map; 130 maps are used spanning the parameter space examined 
here. The results for a convolution kernel 16 pixels wide are shown 
as filled black circles, while open grey circles correspond to a kernel 
with 158 pixels on a side. Different errors for the same magnifica¬ 
tion (vertical spread of the points) are produced due to different 
neighbouring pixels within the range of the convolution kernel, e.g. 
near a caustic. The solid line (red online) is the approximation of 
Equation 0; values of (//) and (N) from any of the maps produce 
practically the same curve. 


smallest profile - the fluctuations disappear and the light 
curve becomes smoother. 

The error in the magnification /j, ;7 of a given map pixel 
is 

, (5) 

where (N) is the average number of rays per map pixel, 
and (//) is the average magnification per pixel, which 
can be computed by the total number of rays in a 
map and the number of rays that would have reached 
the map if there was no lensing (see the Appendix of 
IVernardos fc Flul5 120131 for more details). For con¬ 
volved maps, this error propagates through the convo¬ 
lution formula: 


m,n=—s 

(v) 


_ /]\[\ k m ,n\/Ni +m ,j+n ( 6 ) 

' ' m — —c 


where we have assumed a square, normalized kernel fc, 
with a size of 2s pixels, centered over the i, j-th map 
pixel. Keeping track of the exact error of the convolved 
magnification values means convolving the “error” map 
with the source profile and storing an additional value 
for each light curve pixel. Such an approach would dou¬ 
ble the number of convolutions performed and data that 
have to be stored. To avoid this, we approximate the er¬ 
ror given by Equation ([6]) by assuming that the Poisson 
statistics hold for the convolved ray-count map. In this 
case, the error for the convolved magnification is: 



where ND is not necessarily an integer ray-count any¬ 
more. Taking into account that = N' Lj we 


end up with the final expression for the convolved light 
curve error: 





( 8 ) 


In other words, we approximate the error propagated by 
the convolution with an assumed Poisson error of the 
convolved pixel in question. 

In Figure [31 we compare the above approximation to 
the actual error, calculated in randomly sampled pix¬ 
els from a representative set of maps. We see that our 
approximation is in fact a maximum error of the magni¬ 
fication, which originates from the observation that 
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E 




m,n——s 


E 

m,n=—s 


fcm,n 




( 9 ) 

This inequality depends on the actual convolution ker¬ 
nel used, however, and proving it is beyond the purpose 
of this paper. It suffices to say here that the Gaussian 
profiles we are using are peaked at the location of the 
i, j-th pixel; if profiles that are peaked away from the 
central pixel are used, then our approximation is most 
likely going to fail. 

Thus far we have assumed a fixed physical size for the 
profiles, and the convolutions have been performed for 
a fixed value of I?Ein (see Section However, what 
if a different value of I?Ein is used, e.g. to study spe¬ 
cific lensed systems? In this case, the GERLUMPH light 
curve data can still be used, but caution has to be taken 
on how to scale the source sizes correctly. For example, if 
we apply the GERLUMPH maps to quasar Q2237+0305 
that has Rv,\ n = 1.81 x 10 17 cm for 1 M 0 microlenses 
dMosauera fc KochaneklfeOllIf . our largest source profile, 
1.7 x 10 17 cm, would correspond to a 376-pixel Gaussian 
kernel. The closest convolution kernel to this size is 392 
pixels, which corresponds to a 5 x 10 16 cm profile and 
7?Ein = 5.11 x 10 16 cm. Of course, the closer the size of 
the source profile in pixels is to one of the convolution 
kernels the more similar the light curves will be. The 
case with equal kernel sizes would lead to identical light 
curves, meaning that the scaled light curves would be 
exact, although a different value of i?Ei n was used. 

Variations of the value of i?Ei n due to other reasons, 
e.g. uncertainties in the redshifts of lens and source or 
different mass of microlenses, lead to similar scaling of 
the source profile sizes. We demonstrate this with an ex¬ 
ample of a variation in the Hubble constant. The value 
of i? Ein used in this paper, i.e. 5.11 ± 1.88 x 10 16 cm 
(jMosauera fc Kochanekl[2011f l. was obtained assuming a 
Universe with fl m = 0.3, Q a = 0.7, flk = 0, and Hq = 72 
km s” 1 Mpc -1 . The angular diameter distances appear¬ 
ing in Equation Jl]), Da, depend linearly on the line-of- 
sight comoving distance^ Dc , as 


D a = 


Dc 

1 + z 


( 10 ) 


5 This equation is for an object at redshift z and S i/. = 0. Get- 
ting the angular diameter distance Da-i ,2 between two objects 
at redshifts z\ and Z 2 depends on the geometry of the Universe 
and in general Da-, 1,2 ^ Da -2 ~ Da- i- However, for = 0, 

D A \\,2 = 


1+Z2 
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which in turn is proportional to the Hubble distance, 
Dh = c/Hq. This leads to 

J^Ein _ , - \ 

that relates the error in Hq with the error in the de¬ 
rived i?Ein value. Our profile sizes can then be scaled 
according to the discussion in the previous paragraph. 
Similar reasoning could lead to a different relation for 
the dependence of i?Ein on H 0 for different cosmological 
parameters. 

4. ETOOLS 

4.1. Accessing the data 

The GERLUMPH light curve data can be openly ac¬ 
cessed online. To this purpose, we have extended the 
getquery eTool described in detail in lVernardos fe Flukel 
(I2014al Section 4 and Figure 7). Detailed instructions 
and help tips on how to download the data are provided 
on the GERLUMPH website@ throughout the various 
stages of the process. 

The 2000 light curves for each map are stored in a 
binary file (lc_data.bin), which is compressed using a 
standard Unix tool (e.g. gzip or bzip2). The down¬ 
loaded files are grouped in indexed directories for each 
map and source profile. The values of (/z) and ( N }, re¬ 
quired to calculate the errors in the magnification values, 
are stored in a metadata file (mapmeta.dat) in the map 
directories. Finally, at the root directory there are two 
reference hies for looking up the map and profile indices 
(mINDEX.txt and pINDEX.txt) and their properties, e.g. 
k. 7 , profile size, etc. 

4.2. The GIMLET tool 

Simulated light curves from the GERLUMPH maps 
can be inspected using the GERLUMPH Interactive Mi- 
crolensing Lightcurve Extraction Tool (GIMLET): 

http://gerlumph.swin.edu.au/tools/gimlet/ 

which is openly available online. This is an exploration 
and planning tool, whose main goals are: 

1. To plot high resolution light curves and show their 
location on the corresponding magnification maps. 

2. To plot low resolution light curves for comparison, 
which can be extracted interactively in real time. 

3. To demonstrate the effects of light curve sampling. 

The interactive light curve is not extracted from the full 
resolution magnification map, as this would consume a 
lot of computational resources due to the size of each 
magnification map (381 MB). Instead, a scaled-down ver¬ 
sion of the map is used that has been precomputed and 
stored i n the GER L UMP H datab ase (the sample.png 
file, see lYernardos fe FlukdlWl4all . The spatial resolu¬ 
tion of this map icon is set to 1000 x 1000 pixels, but it 
can be reduced according to the resolution of the monitor 
used to view the webpage. The magnification, in units 

6 visit http://gerlumph.swin.edu.au/guide/ and 
http://gerlumph.swin.edu.au/getquery/ 


of Amag (Equation [3j), is binned in 256 bins in the range 
[—4,4] (values outside this range, which are extreme and 
rare, are set equal to the range limits). These approxi¬ 
mations are used to give a real-time feeling to generating 
light curves. To obtain high quality data users are ad¬ 
vised not to use the GIMLET tool and instead download 
the high resolution light curve data. 

In Figure [4] we show a screenshot of the GIMLET tool 
webpage. The panel shown, “Controls”, contains the ba¬ 
sic features and functionality. High and low resolution 
light curves can be displayed and compared. The length 
of the light curves is measured in units of i?Ein, but it 
is possible to change to units of time (days) by specify¬ 
ing a value for i?Ein and for the effective velocity of the 
source, v s . A time interval, At, can be set to simulate a 
light curve cadence of observations. Other functions can 
be performed in the remaining panels: any of the GER¬ 
LUMPH maps can be loaded (“Change map”), the source 
profile for the interactive low resolution light curves can 
be modified interactively (“Interactive Profile”), and the 
color scheme used on the displayed map can be changed 
(“Color”). 

4.2.1. Implementation 

The content of the tool webpage is rendered by a 
PHP script that handles all the communication with the 
database. JavaScript functions and the jQuerjQ library 
are used to further manipulate the elements of the web¬ 
page: reload the high resolution light curve data, extract 
the interactive light curve, and handle all the button 
functions. The Flot@ library is used to plot the light 
curves. The map is rendered in a HTML50 canvas ele¬ 
ment and the Kineticlibrary is used to plot the light 
curve location on the map. The color is applied to the 
map pixels using the WebGlQ JavaScript application 
programming interface, that allows for GPU acceleration 
at the user end. 

5. SUMMARY AND DISCUSSION 

We have used the GERLUMPH magnification maps to 
produce > 2.5 billion simulated microlensing light curves. 
Our data are publicly available for download from the 
GERLUMPH server. We also release an online explo¬ 
ration and planning tool for plotting the high resolution 
light curves presented here and extracting interactive low 
resolution light curves in real time. Our goal is to provide 
an extensive and consistent set of light curves to be used 
as a benchmark for future parameter space and individ¬ 
ual system microlensing studies. We described numerical 
errors of our data and scaling of the source profiles with 
respect to values of i?Ein and Ho in Section [3] 

Such a complete set of light curves is of high relevance 
to microlensing studies of large numbers of lensed quasars 
in the upcoming all-sky survey era of astronomy. While 
our light curve data cannot be used to fit observed light 
curves for single systems, they are designed to test the 
robustness of, and degeneracies in, such techniques for 
many individual objects in the parameter space. In this 

7 http://jquery.com/ 

8 http://www.flotcharts.org 

9 http://www.w3.org/TR/html5/ 

19 http://kineticjs.com/ 

11 Web Graphics Library: www.khronos.org/webgl 
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Figure 4. Screenshot of the GIMLET tool, located at http://gerlumph.swin.edu.au/tools/gimlet Detailed instructions and help in 
using the GIMLET tool can be found online. Here we outline its main features: (A) Scaled down version of a GERLUMPH map, with the 
locations of the high and low resolution light curves. The low resolution light curve can be moved and rotated interactively by the user. 
(B) The corresponding high (black) and low (grey) resolution light curve variations, plotted as a function of length, measure d in RvAn, or 
time measured in days. (C) Panels that allow the map, source profile, and color coding of the map to be changed (see Section l4.2l for more 
details). (D) The control panel for displaying the light curves and enabling sampling. 


way, we will be able to discover potential systematic er¬ 
rors introduced in the modelling process of microlensed 
quasars, that will hopefully lead to better measurements 
of Hq and accretion disc constraints. 

Our results can be used to train ma chine learn i ng algo¬ 
rithms for measu ring time delays (e.g. lCuevas-Tello et al.l 
120061 iHirv et all 120111) . Our simulations do not try to 
reproduce any specific observed light curve directly - in¬ 
stead, they cover a large range of possible, yet currently 
unobserved scenarios. As such, they are highly suitable 
as an input to unsupervised machine learning algorithms 
as part of the process of automatically discovering and 
classifying thousands of new lensed quasars set to be dis¬ 
covered in future synoptic surveys. For example, they 
are ideally suited for the time de lay challenge (TDC; 
iDobler et al.ll2013i iLiao et ahll2014D , a collaborative com¬ 
munity effort that uses mock observations of multiply im¬ 
aged quasars and attempts to measure the time delays 
using various techniques. The GERLUMPH light curves 
are available at high resolution and for a wide range of 
source and lens model parameters that can be matched 
to those of the TDC mock observations. 

Additional machine learning and data mining ap¬ 
proaches can be use d to explore our simulated light 
curve dataset (e.g. see lBall & Brunncrll2010l : Hvezic et alJ 
120141 ). This could be done by calculating basic statisti¬ 
cal properties, e.g. mean, median, variation, etc, or more 
advanced properties, like the power spectrum, the num¬ 
ber of peaks in the light curves, their prominence and/or 
separation, etc. The entire sample of light curves could 
then be classified using such a metric as a new way of 


exploring and understanding lens and source model spe¬ 
cific degeneracies. Given the large size of the data, i.e. 
2.5 billion light curves corresponding to 5.6 TB of data, 
the high dimensionality of the parameter space, and the 
many metrics and classification techniques that can be 
used, this task will be the topic of future work. 

We have focused on the size as the most significan t 
source factor for microlensing (jMortonson et al1l2005l ). 
However, theoretical studies could be envisaged using 
more complicated source profiles (inclined discs, bicon- 
ical flows, disc with hot spots, etc) that could be com¬ 
pared against the dataset presented here. Since our fast 
GPU convolution implementation (> 1.28 million convo¬ 
lutions across 5 days using gSTAR) and our data man¬ 
agement infrastructure are already in place, this is a rel¬ 
atively straightforward task. Further developments, like 
machine learning classification techniques, would provide 
additional tools to investigate the effect of second order 
source characteristics. 

Our choice of 1.5 i?Ein~long light curves, and 2000 light 
curves per source and map was justified in Sections [2] 
and [3] respectively. The main restriction on the length 
and number of light curves is the storage space currently 
available on the host facility. Whereas much of the mi¬ 
crolensing analysis process (to date) has been compute 
or memory limited, we are now in a regime where we are 
I/O limited. We are separately investigating the useful¬ 
ness of data compression techniques (lossless and lossy) 
which ma y allow us to add additional light curves to the 
database (lYohl et al.ll2015l) . 

Planning of future and ongoing monitoring campaigns 
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of specific systems could be facilitated by using our high 
resolution data and the GIMLET tool presented in Sec¬ 
tion [4] Existing GERLUMPH maps for the k, 7 values 
of the targeted systems can be selected from the GER¬ 
LUMPH database and inspected. By providing a value 
for i?Ein and the effective source velocity, low resolution 
light curves can be produced interactively by changing 
their length and orientation on the map. The user can 
then intuitively decide the best observational strategy, 
e.g. in terms of the cadence of observations, based on 
the visual appearance of the light curves. 

The GIMLET tool presented in this paper is a building 
block for a more advanced online modelling tool. GPU 
acceleration in the browser provided by WebGL presents 
opportunities to model three-dimensional source profiles 
interactively in real time, e.g. to add components like 
an event horizon or jets, and modify their properties, 
like the opening angle or Schwarzschild radius. WebGL 
could also be used to perform the convolutions between 
the modelled source profiles and the scaled down version 
of the map without any overhead, since all the computa¬ 
tions would be performed in the browser in real time. Fi¬ 
nally, we intend to integrate the browser-based front-end 
more fully with a back end supercomputer, where user 
requests for high resolution full scale modelling would be 
sent for computation. This working model may in fact 
not lie far from data access and analysis approaches that 
will be followed by the future all-sky survey facilities. 

In conclusion, we have presented high quality mi- 
crolensing simulations throughout the parameter space, 
preparing the theoretical ground for the upcoming all 
sky survey era of quasar microlensing. This dataset can 
be used as a benchmark for existing and future single 
object and parameter space studies, and can be further 
explored using machine learning and data mining tech¬ 
niques. Our data and software is made openly available 
for further use by the community. We complement our 
data with comprehensive and innovative online analysis 
tools. We hope that this work will contribute to the ad¬ 
vancement of data intensive quasar microlensing studies 
of the future. 

This research was undertaken with the assistance of 
resources provided at the GPU-Supercomputer for The¬ 
oretical Astrophysics Research (gSTAR) through the As¬ 
tronomy Supercomputer Time Allocation Committee, 
supported by the Australian Government. gSTAR is 
funded by Swinburne University and the Australian Gov¬ 
ernment’s Education Investment Fund. DC acknowl¬ 
edges receipt of a QEII fellowship from the Australian 
Research Council. 
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